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Spiral shocks in the accretion disc of IP Peg during 
outburst maximum 
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ABSTRACT 

In response to our recent discovery of spiral arms in the accretion disc of IP 
Peg during rise to outburst, we have obtained time-resolved spectrophotometry of 
IP Peg during outburst maximum. In particular, indirect imaging of He II 4686, using 
Doppler tomography, shows a two-arm spiral pattern on the disc image, which confirms 
repeatability over different outbursts. The jump in He II intensity (a factor of more 
than two) and in velocity (~200-300 km s _1 ) clarifies the shock nature of the spiral 
structure. The He II shocks show an azimuthal extent of ~90 degrees, a shallow power- 
law emissivity ~ an upper limit of 30 degrees in opening angle, and a flux 
contribution of 15 per cent of the total disc emission. We discuss the results in view 
of recent simulations of accretion discs which show that spiral shocks can be raised in 
the accretion disc by the secondary star. 
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1 INTRODUCTION 

High-energy phenomena such as jets emanating from AGNs 
and X-ray radiation emitted from binaries are related to 
accretion discs. Accretion discs are the most efficient ma- 
chines for the extraction of energy and angular momentum. 
The physical conditions that lead to specific dynamical for- 
mations in discs, such as spiral arms, are not contained as 
yet due to lack of observational constraints. Spiral structure 
in galaxies is attributed to a quasi-stationery wave which 
triggers star formation as it propagates through the disc 
(Bertin et al. 1989) or a tidal pattern due to interaction 
with a satellite galaxy (e.g. M51; Toomre and Toomre 1972). 
Spiral waves in a disc have also been invoked to explain the 
proximity of giant extra-planets to their orbiting stars (e.g 
51 Peg; Lin et al. 1996). Moreover, spiral shocks have been 
found in simulations of protosolar discs and have been pro- 
posed as a means of converting the gaseous protostellar disc 
into orbiting planetesimals (Boss 1997). Accretion discs in 
cataclysmic variables (a white dwarf accretes matter from a 
donor K-M star) evolve on short dynamical timescales (2-10 
hours) allowing us to investigate in detail their development. 
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In cataclysmic variables, the outburst origin in dwarf novae 
(disc or donor star) and the mechanism for the angular mo- 
mentum transport of the disc material ('viscosity') are still 
unresolved issues though fundamental in our understanding 
of accretion physics (Verbunt 1986). The exact structure of 
the disc, an a-disc (Shakura and Synyaev 1973) or a spiral- 
shock disc (Sawada, Matsuda & Hachisu 1986), is linked to 
the above issues. Simulations of such accretion discs have 
revealed double spiral shocks (Sawada, Matsuda & Hachisu 
1986; Savonjie, Papaloizou and Lin 1994; Stehle 1998). Re- 
cently, we indirectly detected spiral waves in the accretion 
disc of the eclipsing dwarf nova IP Peg during the rise to 
outburst (Steeghs, Harlaftis, Home 1997). Here, we report 
on subsequent observations of He II 4686, during maximum 
of the November 1996 outburst, which were aimed to probe 
the ionization structure of the spiral arms. 



2 OBSERVATIONS 

IP Peg, a double-eclipsing dwarf nova, shows semi-periodic 
outbursts every ^3 months. During the third day of the 
November 1996 outburst (Fig. 1), we observed a complete 
binary cycle of 3.8 hours. We obtained 81 spectra, under 
« 1.6 arcsecond seeing, with a TEK CCD and the 235 
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Figure 2. Doppler tomography of He II 4686. From left to right, the panels display the observed data, the Doppler image, and the 
projection of the Doppler image into a trailed spectrogram. The image shows non-axisymmetric components arising from the red dwarf, 
the disc (tidal shocks) and an unresolved location with almost zero velocity. The red dwarf size and ballistic trajectory are also drawn 
(for K c = 280 and K w = 162 km s _1 ). Steps in units of Rl 1 are marked along the ballistic velocities. The greyscale bar corresponds to 
a linear intensity scale of 0-57 mJy. 
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Figure 1. The optical light curve, as given by AAVSO, during 
the November 1996 outburst of IP Peg. Eclipses are visible as 
sharp drops in magnitude. The time of the INT observations is 
marked. 



and sky-subtraction. The wavelength calibration was per- 
formed using 17 CuAr arc lines and a 3rd order polyno- 
mial fit and is accurate to 0.04 A. A comparison star had 
been included in the slit which was subsequently used to 
correct the object spectra for atmospheric and slit losses. 
The absolute flux scale was denned by observing the flux 
standard Feige 15 (Oke 1990). Typical signal to noise ra- 
tio is 21 for the IP Peg spectra (and 46 for the compari- 
son spectra). Flux uncertainties were estimated using Pois- 
son statistics and a noise model for the CCD chip (gain of 
1.1 electrons per ADU and readout noise of 5.7 electrons 
per pixel). We adopt the binary ephemeris from Wolf et 
al. (1993) T (HJD) = 2445615.4156(4) + 0.15820616(4) E, 
where T Q is the inferior conjunction of the white dwarf. The 
He II A4542/A4686 flux ratio of 0.022-0.036, as well as the 
He I flux ratio A6678/A4471 flux ratio of 0.71-0.87 (see also 
Marsh and Home 1990), are consistent with recombination 
case B (for T = 5,000-20,000 K and electron density of 10 4 
cm -3 ; Osterbrock 1989). Given that the He I is blended with 
the Mg II line, the above ratio is only a lower limit, thus 
favouring the temperatures nearer to 20,000 K. 



camera of the Intermediate Dispersion Spectrograph on the 
2.5m Isaac Newton Telescope at La Palma. The wavelength 
range covered is 4354-4747 A at a velocity dispersion of 
27 km s _1 per pixel. The spectra were extracted using op- 
timal extraction (Home 1986) after debiasing, flat-fielding 



3 RESOLVING SPIRAL WAVES WITH 
DOPPLER TOMOGRAPHY 

The trailed spectra of the strongest emission line, He II, show 
very complex structure (left panel in Fig. 2). The peak sepa- 
ration changes with phase and is narrower at phase 0.50. The 
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Figure 3. A model using spiral arms (top-right) is shown for comparison to the observed non-axisymmetric image (top-left). The 
computed data resulting from the model projected at each binary phase (bottom-right) are also displayed for comparison with the 
non-axisymmetric part of the data (bottom-left). The signature of spiral shocks is demonstrated with clarity on the trailed spectra (a 
decreasing double-peak separation with an intensity jump at phase 0.6). The greyscale bar corresponds to a linear intensity scale of 0-100 
for the images and of 0-25 mJy for the data. 



red peak of the profile is maximum at phases 0.45 and 0.95 
where it has a velocity of 400 and 650 km s _1 , respectively. 
Similarly, the blue peak has a clear maximum at phase 0.85 
and a velocity of 500 km s _1 . Other weaker components are 
immediately visible; a sharp component moves from red to 
blue at phase 0.5 (red star emission) and a low- velocity com- 
ponent is present at phases around 0.2 and 0.9. We recon- 
struct the Doppler image of the He II 4686 emission-line dis- 
tribution using the trailed spectra (Marsh and Home 1988). 
This imaging technique has been particularly successful in 
resolving the location of emission components such as the 
red star (IP Peg; Harlaftis et al. 1994), the gas stream (OY 
Car in outburst; Harlaftis and Marsh 1996), the bright spot 
(GS2000+25; Harlaftis et al. 1996) and spiral waves in the 
outer accretion disc (from Ha and He I lines of IP Peg during 
rise to outburst; Steeghs, Harlaftis and Home 1997). 

The He II Doppler image of IP Peg during outburst 



(centre panel of Fig. 2) directly displays the locations of 
the various emission components; the inner side of the red 
star, a low-velocity component and the dominant accretion 
disc with extended spiral arms. For comparison with the 
observed data, the computed data from the image are also 
shown (right panel). The low- velocity component is centred 
at (V x ,Vy) = (-100 ± 50,-20 ± 70) km s _1 and has a 
FWHM of 270 km s _1 and although is seen in other dwarf 
novae during outburst its origin is not understood (Steeghs 
et al. 1996). In Fig. 3 we display the non-axisymmetric data 
and Doppler image for comparison with a model. The non- 
axisymmetric Doppler image (top-left) is built by subtract- 
ing the median at each radius with the white dwarf at its 
centre. A simple model of spiral arms (plus a red star and a 
low- velocity component; top-right) is built by using a power 
law of V~ x for the intensity, velocities between 400 to 800 
km s _1 and azimuthal extent of 0.25 binary cycles. The com- 
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Figure 4. The variation of maximum intensity and velocity at 
maximum intensity with binary phase, as extracted from the He II 
Doppler image (Fig. 2) (see also text). 

puted data from the model (bottom-right; with the same 
noise content as the data) reconstruct successfully all the 
velocity features observed, such as the decreasing double- 
peak separation with an intensity jump at phase 0.6. The 
Bowen blend, He I 4471 and Mg II 4481 lines produce simi- 
lar Doppler images. 



4 DISCUSSION 

The high-ionization He II 4686 line is most suitable to clar- 
ify the shock nature of the resolved spiral structure, since it 
is sensitive to higher temperatures and is not saturated as 
Ha may be. We quantify the properties of the spiral struc- 
ture by taking radial slices starting from the white dwarf at 
(0,-147) km s _1 for each azimuth and then search for the 
velocity for which the image intensity reaches its maximum 
value. The result is then to trace the centre of the spiral 
arms with azimuth, maximum intensity (bottom) and ve- 
locity at that image value, relative to the white dwarf (top). 
The two peaks in both velocity and intensity correspond to 
the spiral arms. In particular, the jump in intensity (a fac- 
tor of more than two) and velocity (^200-300 km s _1 ) is 
the characteristic signature for a shock. The shocks show 
an azimuthal extent of ~90 degrees (phases 0.50-0.75 and 
0.95-1.25) and a very shallow power- law index in emissivity 
(rsj V~ x or rsj R~ 2 5 assuming a Keplerian disc). The spiral 
shocks contribute about 15 per cent of the total disc emis- 
sion and their opening angle is ~ 30 degrees by assuming 
Keplerian flow (which places a firm upper limit). The Mach 
numbers at the shocks range between 28-50 for T — 20, 000 
K and velocities between 400-700 km s -1 . There are differ- 
ences between the two arms: the arm centred at phase 0.65 
(bottom-left arm on Doppler image) is weaker in intensity 
(by 40 per cent) than the arm centred at phase 0.05. Also, 
the velocity range is 400-700 km s" 1 compared to 400-550 
km s _1 and the radial width is 40 per cent shorter compared 
to the FWHM=240 km s" 1 of the arm at phase 0.05. Note 
that in the above we have discussed the velocities related to 
the maximum intensities of the image. The Doppler image 



shows that there is also some emission at velocities up to 
1000 km s _1 (or 850 km s _1 relative to the white dwarf). 

Spiral arms have now been observed in 3 different 
outbursts of IP Peg (this work; 1993 August outburst in 
Steeghs, Harlaftis and Home 1997; 1987 July outburst in 
Marsh & Home 1990 where, in retrospect, only limited ev- 
idence for spiral structure can be discerned due to the in- 
complete binary coverage). The trailed spectra obtained by 
Hessman (1989) immediately after the end of an outburst of 
IP Peg do not show any structure other than the red star 
'S'-wave component. We tested this by adjusting our He II 
spectra to the instrument and time resolution of the Hess- 
man spectra (195 km s _1 and 12 phase bins, respectively). 
The signal-to-noise is similar since the red star line compo- 
nent is visible. The spiral shocks can still be resolved in the 
simulated trailed spectra suggesting that the shocks either 
are tightly wound up or may have dissipated by the time 
the system has reached quiescence after an outburst. Spiral 
structure in systems other than IP Peg can also be found, 
such as in SS Cyg during outburst (see Helium lines in Fig- 
ures 7, 9 and 10; Steeghs et al. 1996) and in the Ha trailed 
spectra of LB1800 (Still et al. 1998, Steeghs et al. 1998, 
in preparation). No evidence of spiral shocks is observed in 
short-period SU UMa-type stars, such as OY Car during 
outburst (Harlaftis and Marsh 1996), SU UMa and YZ Cnc 
in outburst or quiescence (Harlaftis 1991). Therefore, it may 
be possible that the spiral shocks only develop in systems 
above the period gap where the disc, the companion star and 
the mass transfer rate are larger. Doppler maps of X-ray bi- 
naries do not show any evidence of spiral shocks (outburst 
of GRO J0422+32, Casares et al. 1995; quiescent map of 
A0620-00, Marsh et al. 1994; quiescent map of GS 2000, 
Harlaftis et al. 1996). Only Ra spectra of X 1822-371, the 
prototype system for vertical disc structure, may contain a 
hint of some structure similar to spiral shocks (Harlaftis et 
al. 1997). 

SPH simulations of the IP Peg accretion disc do show 
the development of transient spiral shocks (Armitage and 
Murray 1998) which agree quite well with the observed ones 
in velocity and azimuthal extent. Stehle (1998) uses hy- 
drodynamical simulations of a hot accretion disc (T e // = 
10, 000 - 50, 000 K,q = 0.3) which confirm the deployment 
of stable spiral shocks in the disc over many orbital cycles 
(see also Steeghs and Stehle 1999). Fitting of the He II spiral 
shock data (Fig. 4) with model data extracted from simu- 
lated discs should eventually constrain the physical condi- 
tions of the disc such as the local Mach number and effective 
temperatures. We are aiming to obtain similar coverage at 
the end of an outburst (to complete our coverage), so that 
we have a full picture of the dynamical evolution of the tidal 
shocks. In particular, the effect of a decreasing a parameter 
(outburst decline) on the velocity range, opening angle and 
azimuthal extent of the spiral shocks. 
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